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We study the stability margin of a vehicular formation with distributed control, in which the control 
at each vehicle only depends on the information from its neighbors in an information graph. We consider 
a D-dimensional lattice as information graph, of which the 1-D platoon is a special case. The stability 
margin is measured by the real part of the least stable eigenvalue of the closed-loop state matrix, which 
quantifies the rate of decay of initial errors. In [1], it was shown that with symmetric control, in which 
two neighbors put equal weight on information received from each other, the stability margin of a 1-D 
vehicular platoon decays to as 0(1/ N 2 ), where N is the number of vehicles. Moreover, a perturbation 
J> 1 analysis was used to show that with vanishingly small amount of asymmetry in the control gains, the 

stability margin scaling can be improved to 0(1 /N). In this paper, we show that, with judicious choice 
of non-vanishing asymmetry in control, the stability margin of the closed loop can be bounded away 
qq ■ from zero uniformly in N. Asymmetry in control gains thus makes the control architecture highly 

scalable. The results are also generalized to D-dimensional lattice information graphs that were studied 
in 0, and the correspondingly stronger conclusions than those derived in J2] are obtained. In addition, 
we show that the size-independent stability margin can be achieved with relative position and relative 
velocity (RPRV) feedback as well as relative position and absolute velocity (RPAV) feedback, while the 
analysis in Q], (2) was only for the RPAV case. 
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I. Introduction 

We study cooperative control of a large vehicular formation with distributed control. The 
vehicles are modeled as double integrators, and the control action at each vehicle is computed 
based on information from its neighbors, where the neighbor relationship is characterized by 
a lattice information graph. The control objective is to make the vehicular formation track 
a constant-velocity type desired trajectory while maintaining pre-specified constant separation 
among neighbors. The desired trajectory of the entire vehicular formation is given in terms of 
trajectories of a set of fictitious reference vehicles. 

The problem of distributed control for multi-agent coordination is relevant to many applications 
such as automated highway system, collective behavior of bird flocks and animal swarms, 
and formation flying of unmanned aerial and ground vehicles for surveillance, reconnaissance 
and rescue, etc. fl3j — [[81 - A typical issue faced in distributed control is that as the number of 
agents increases, the performance (stability margin and sensitivity to external disturbances) of 
the closed loop degrades. Several recent papers have studied the scaling of performance of 
vehicle formations as a function of the number of vehicles. The references 0J, have studied 
the scaling of the stability margin of D-dimensional lattice formations. The stability margin is 
defined as the absolute value of the real part of the least stable eigenvalue of the closed loop. The 
stability margin characterizes the rate at which initial errors decay. The references ll9ll- [fT3l have 
examined the sensitivity of 1 -dimensional platoons to external disturbances. However, among 
papers that examined sensitivity to disturbance, to the best of our knowledge only [13J has 
considered asymmetric control, the rest are limited to symmetric control. The control is called 
symmetric if between two neighboring vehicles i and j, the weight i puts on the information 
from j is the same as the weight j puts on the information from i. 

In previous works on 1-D vehicular platoons, two types of feedback are respectively con- 
sidered: relative position absolute velocity (RPAV) feedback |[T], lfl2l and relative position 
relative velocity (RPRV) feedback |fTT|. 0"3ll . 031. With symmetric control, the stability margin 
of the vehicular platoon decays to as 0(1/ TV 2 ) in both types of feedback. This result for 
RPAV feedback was shown in [1], and for RPRV feedback was shown in [fl4l . The loss of 
stability margin with symmetric control has also been recognized by other researchers lfl2ll . 
031. Asymmetric control in the RPAV case was examined in 0}, 0, where it was also shown 
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that with vanishingly small asymmetry in the control gains, the stability margin can be improved 
to 0(1/N). Similar conclusions are also obtained for a vehicle formation with a D-dimensional 
lattice as its information graph - that decay of stability margin can be improved with 
asymmetry. In case of RPRV feedback, a similar improvement to 0(1/ N) with asymmetry 
was shown in |fl4l . where only the relative velocity feedback gains were made asymmetric. 
The analyses in flU, 0, 031 were based on a partial differential equation (PDE) approximation 
of the closed loop dynamics and a perturbation method; the latter limited the results to only 
vanishingly small asymmetry. 

In this paper we provide a stronger result on the stability margin with asymmetric control 
by avoiding the perturbation analysis of the aforementioned papers. We also avoid the PDE 
approximation and analyze the state space model directly. In particular, we show that with 
judicious choice of asymmetry in the control, the stability margin of the vehicular formation can 
be uniformly bounded away from (independent of N) and derive a closed-form formula for 
the lower bound. This result makes it possible to design the control gains so that the stability 
margin of the system satisfies a pre-specified value irrespective of how many vehicles are in the 
formation. We also generalize the result to formations with D-dimensional information graphs, 
and show that a similar, size-independent stability margin can be obtained by using asymmetry 
in the control gains. These results are established for both RPAV and RPRV feedbacks. 

The focus of this paper is on the stability margin, which is related to exponential stability 
of the closed loop system. A related concept is that of "string stability" lfT6ll . String stability is 
usually interpreted as the system's sensitivity to external disturbances; see J6]|, ifTOll . ifTTll . lfT8l 
and references therein. We do not study sensitivity to external disturbances in this paper. 

For ease of description, we first present the problem statement and main result for a vehic- 
ular formation with 1 -dimensional information graph (i.e. the vehicular platoon) in Section UH 
Analysis of the stability margin and numerical verification appear in Section Hill The extension 
of the result to a vehicular formation with D-dimensional lattice information graph is presented 
in Section [IV] The paper ends with a summary in Section |V] 
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Fig. 1. Desired geometry of a vehicular platoon with N vehicles and 1 "fictitious" reference vehicle. The filled vehicle in the 
front of the platoon represents the reference vehicle, it is denoted by index 0. 



II. Problem statement and result for 1-D platoon 
A. Problem statement 

In this section we consider the formation control of N homogeneous vehicles which are moving 
in 1-D Euclidean space, as shown in Figure \T\ The position of the i-th vehicle is denoted by 
Pi E K and the dynamics of each vehicle are modeled as a double integrator: 

pi = Ul , i E {1, 2, • • • , N}, (1) 

where Ui E K is the control input. This is a commonly used model for vehicle dynamics in 
studying vehicular formations, which results from feedback linearization of non-linear vehicle 
dynamics Ifl9i Il20l. 

The control objective is that vehicles maintain a desired formation geometry while following 
a constant-velocity type desired trajectory. The desired geometry of the formation is specified by 
the desired gaps A^-i^) for % E {1, • • • , N}, where A^-i^) is the desired value of Pi-i(t) —pi(t). 
The desired inter- vehicular gaps A^-i^'s are positive constants and they have to be specified in 
a mutually consistent fashion, i.e. A( i)fe ) = A^+Ay ^ for every triple (i, j, k) where % < j < k. 
The desired trajectory of the platoon is provided in terms of a fictitious reference vehicle with 
index 0, whose trajectory is given by Po(t) = v*t + Co for some constants v*,co, where v* is 
the cruise velocity of the formation. The desired trajectory of the i-th vehicle, p*(t), is given by 

?m = pin) - a (w) = P * (t) - e; =1 a w _ w) . 

We consider the following distributed control laws. 

1) Relative position and absolute velocity (RPAV) feedback: the control action at the i-th 
vehicle depends on the relative position measurements with its two neighbors (one on either 
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side), its own velocity, and the desired velocity v*: 

Ui = - k{(pi - + A (l _i >i) ) - k\{pi - p i+ i - A( iji+1 )) 
-biipi-v*), i E {1, ■ ■ ■ , N — 1}, 

UN = ~ k f N (p N ~ PN-1 + A(jV-l,jV)) - ~ ( 2 ) 

where A;/, A;* are the front and back position gains and hi is the velocity gain. 

2) Relative position and relative velocity (RPRV) feedback: the control action at the i-th vehicle 
depends on the relative position and relative velocity measurements with its nearest neighbors 
in the platoon: 

Ui = - k{(pi - p^i + A (i _i ii} ) - k\{pi - p. i+ i - A (M+ i)) 

- b{(pi - pi-!) - - p i+1 ), i E {1, ■ • • , N - 1}, 

u N = - k N (p N - Pn-i + A(jv-i,jv)) - b f N (p N - PN-l), (3) 

where k{, k\ (respectively, b{, b\) are the front and back position (respectively, velocity) gains 
of the i-th vehicle. 

In the RPRV feedback case, vehicle i must be provided (a-priori) the desired gaps with its 
two neighbors. In the RPAV feedback, it must be provided with additional information: the 
formation's desired velocity v*. The closed-loop dynamics with RPAV (resp., RPRV) feedback, 
in terms of the tracking errors pi '■ = Pi — p*, can now be expressed as: 

x = A {RFAV) x, (resp.) x = A (rprv) x, (4) 

where the state vector is defined as x := \p!,pi, ■ ■ ■ ,Pn,Pn] E M? N , and the state matrix 
depends on the control gains but not on the desired gaps or desired velocity. 

Definition 1: The stability margin S^^' (respectively, 5'^ RPRV ^) of the closed-loop system © is 
defined as the absolute value of the real part of the least stable eigenvalue of A( WA!f > (respectively, 
^(rprv)^ rp^ £ contro i i aw Q (respectively, ©) is symmetric if each vehicle uses the same front 
and back control gains: k{ = k\ = k , \ = b (respectively, k{ = k\ = k , b{ = b\ = b ), for all 



i E {1, 2, • • • , N — 1}, where k , b are positive constants. □ 
In this paper, we consider the following asymmetric control gains 

RPAV feedback: k{ = (1 + e)k , k\ = (1 - e)Jfe 0> h = b . (5) 
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RPRV feedback: 



k{ = (l + e)k , k\ 
b{ = (l + e)b , b\ 



(l-e)^o, 



(1 - e)6o, 



(6) 



where e E [0, 1) denotes the amount of asymmetry; e = corresponds to symmetric control. 
The design for the RPAV case is inspired by |0Q, 0. The control gains given in (0) and © are 
homogeneous in the sense that they do not vary with i. The reason we only consider homogeneous 
control gains is that heterogeneity has little effect on the scaling of stability margin, see lfT4l 
for a proof for 1-D platoon. The proof for vehicular formation with general graphs is given in 
Lemma [2l which is provided in the appendix. 

The following proposition summaries the results in flU, Ifl4ll . 

Proposition 1: Consider an iV- vehicle platoon with closed loop dynamics ©. 

1) [Corollary 1 of 0], Theorem 1 of El] With symmetric control (e = 0), both and 



2) [Corollary 3 of HI] With the asymmetric control gains k{ = k (l + e), k\ = k (l — e) and 
bi = b , the stability margin of the platoon with RPAV feedback is 5 (rpav) = O(jj). [ 

3) [Theorem 2 of Ifl4ll 1 With asymmetric control gains k{ = k\ = k ,b{ = b (l + e),b\ = 
b (l - e), the stability margin of the platoon with RPRV feedback is 5"( RPRV ) = O(jj). 



Proposition \T\ shows that with symmetric control, the stability margin decays to as 0{1/N 2 ), 
irrespective of the type of feedback we used. However, in the case of RPAV feedback, with 
vanishingly small amount of asymmetry in the position gains, the stability margin of the system 
can be improved to 0(1/N). The same 0(1 /N) trend can be achieved for the case of RPRV 
feedback with vanishingly small asymmetry in the velocity gains alone while the position gains 
are held symmetric. The design © was not considered in [|T4|. Since the results in JTJ, [[T4| 
were obtained with a perturbation analysis, these results are applicable only when the amount 
of asymmetry is vanishingly small. 

The following theorem is the main result of this paper, whose proof and numerical corrobo- 
ration are given in Section UTT1 

'The case considered in |1J was that \k{ — ko\ < e , |fej — ko\ < e. It is straightforward, however, to re-derive the results if 
the constraints on the gains are changed to the form used here: \k( — ko\/ka < e , \k% — ko\/ko < e. In this paper we consider 
the latter case since it makes the analysis cleaner without changing the results of [1| significantly. 



5 (rprv) ^ ( i ). 



Statements (2) and (3) hold in the limit e — > and N — > oo. 



□ 
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Theorem 1: With the control gains given in © and © respectively, for any fixed e G (0, 1), 
the closed loop is exponentially stable and the stability margin of the vehicular platoon is bounded 
away from uniformly in N. Specifically, 



5 (rpav) >^ y / (7) 



(8) 



where 3ff(.) denotes the real part. □ 
Remark 1: Comparing Theorem \T\ with Proposition [T} we observe the following: (1) Even with 
an arbitrarily small (but fixed and non-vanishing) amount of asymmetry in the control gains, the 
stability margin of the system can be bounded away from zero uniformly in N. This asymmetric 
design therefore makes the resulting control law highly scalable; it eliminates the degradation of 
stability margin with increasing N. (2) In case of the RPAV feedback, although the control law 
is the same as that analyzed in [[TJ, the stronger conclusion we obtained - compared to that in 0J 
- is due to the fact that our analysis does not rely on a perturbation-based technique that was 
used (TJ, which limited the analysis in [[TJ to vanishingly small e. (3) For the RPRV feedback 
case, the stronger result compared to that in 031, is obtained by putting equal asymmetry in 
both position and velocity gains, while [fT4ll allowed asymmetry only in the velocity gain. In 
addition, unlike fl], |fl4ll . we do not use a PDE (partial differential equation) approximation to 
analyze the stability margin, but analyze the state- space model directly. □ 

III. Stability margin of the 1-D vehicular platoon 

With the control gains specified in © and © respectively, it can be shown that the state 
matrices can be expressed in the following forms, 

^(RPAV) = J n ^ ^ + L (l) ^ 

^(RPRV) = j N q As + L (l) ^ ^ (9) 

where I N is the N x N identity matrix, <g> denotes the Kronecker product, and 



1 












, A 2 := 






-6q 




-k 
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^3 := 








A A := 
















-k -b 



(10) 



where k > 0, bo > are the nominal position and velocity gains respectively, and 

2 -1 + e 
-l-e 2 -1 + e 



-1 + e 
1 + e 



(11) 



-l-e 2 
-1 - 

It follows from Theorem 3.1 of [|2TT| that the eigenvalues of are given by 

A = o + 2cpcosfl, (12) 

if 9 (9 7^ m7r, meZ, Z being the set of integers) is a solution to 

p N (acsm(N + 1)5 + (^5 - a/3) sin(iV - 1)9 
-cp(7 + 5) sin NO) - (cap 2N + a/3) sin 9 = 0, (13) 

where a = -l-e,o = 2,c=-l + e,a = /3 = 7 = 0,o" = -l + e,p = - e)/(-l + e). 

Eq. (fT2l) and (TTSl can now be simplified to 

X e = 2-2Vl-e 2 cos0 € , £ G {1,2,--- ,iV}, (14) 

where e G (0, 1) and 9 t is the £-th root of the following equation 



1 



sin(iV + 1)9 = sin N9. 



1 - e 

From (fl4"l) . we see that the eigenvalues of are real and positive, and moreover, < Ai 

3tt 



(15) 



2 - 2\/l - e 2 cos 9i < A 2 < ■ • • < A^ = 2 - 2\/l - e 2 cos0 N , where 0i G 



#/vG 



^2(JV+1) ' 2(JV+1)/' 

( 2(N+l) ; ^(AT+i) ) are me so l ut i° ns t0 (TT3T ) . To see why, first notice that we only need consider 
the roots of (fl"5T) in the open interval (0, 2ir), in which there are 2iV nontrivial isolated roots. See 
Figure |2] for an example. The roots located in R \ (0, 2%) are 2mix (m G Z) distance away from 
those in (0, 2tt). Moreover, if 9 G (0, 2n) is a solution of (fl"5l) . then 2n — 9 is also a solution. 
Therefore, we can restrict the domain of analysis to (0, it), in which there are N isolated roots. 
The ordering of the eigenvalues follows from cos 9 being a decreasing function in (0, 7r). It is 
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Fig. 2. Graphical solution 9 of ^/(l + e)/(l - e) sin((JV + 1)0) = sin(iV0) with e = 0.1 and N = 3. 



straightforward to show from graphical solution of (fl"5T) that the £-th root 6>£ is in the open interval 
( i 2W+ip ' ^ e now P resent a formula for the stability margin of the vehicular platoon in 
terms of the eigenvalues of LP-\ 

Lemma 1: With the control gains given in © and © respectively, and < e < 1, the stability 
margin of the vehicular platoon is 



^i(rpav) 



g(RPRV) 



bo 
2 ' 

bo-y/bl-4k \ 1 
2 



6qAi 
2 ' 

2kp 



if \i > b 2 /4k , 
otherwise, 



if A,v < 4k /h 



o- 



mm 



jjnAi 



2fcp 



if Aa > Ako/b 2 Q , 
otherwise, 



2 ' b + y/bl-Ako/X f 

where Ai and A^ are the smallest and largest eigenvalues of respectively. 



□ 



Proof of Lemma [7] Our proof follows a similar line of attack as of ll22l . From Schur's triangu- 
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larization theorem, there exists an unitary matrix U such that 

U- l L^U = L u , 

where L u is an upper-triangular matrix whose diagonal entries are the eigenvalues A^ of 
We first consider the RPAV feedback case. We do a similarity transformation on matrix A^ AV \ 

= (f/" 1 ® 7 2 )(/ iV ® A + ® A 2 )([/ ® 7 2 ) 
= iiv (8) Ai + L u ® A 2 . 

It is a block upper-triangular matrix, and the block on each diagonal is A\ + X^A 2 , where A<? G 
^•(Z/ 1 )), and cr(-) denotes the spectrum (the set of eigenvalues). Since similarity transformation 
preserves eigenvalues, and the eigenvalues of a block upper- triangular matrix are the union of 
eigenvalues of each block on the diagonal, we have 

a(A^) = a(A^) = |J {a(A 1 + X e A 2 )} 



a 



1 
—k \i —bo 



(16) 



= u 

It follows now that the eigenvalues of A^^ are the roots of the characteristic equation s 2 + 
b s + k X e = 0. For each I G {1, 2, • • • , iV}, the two roots are 

, _ -jg ± y/feg - 4k X e 
b e — 2 ^ ' 

The root closer to the imaginary axis is denoted by sj~, and is called the less stable eigenvalue 

between the two. The least stable eigenvalue is the one closet to the imaginary axis among them, 

it is denoted by s min . It follows from Definition Q] that S = |9ft(s min )|. 

Depending on the discriminant in (TT71 ), there are two cases to analyze: 

1) If Ai > &q/4&o, due to Ai < • • • < A at, we have the discriminant in (fTTT ) for each I is 



non-positive, which yields 5 ,( ^ RPAV - ) = |3?(s r 



6o 
2 • 



2) Otherwise, the less stable eigenvalues can be written as s\ = |(— b + a/^o — Ak^Xi), 
which may be complex for some t > 1. The least stable eigenvalue is obtained by setting 

X £ = X u so that = |»( Smin )| = |(&o - ^o-^oAi)- 
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For the case of RPRV feedback, following the same procedure as that of RPAV feedback, the 
characteristic equations are given by 

s 2 + X £ b s + Xik = 0. (18) 

For each £ £ {1, 2, • • • , N}, the two roots of the characteristic equations (fT8l) are, 

± _ Xjbp _^ a/ (X e b ) 2 - 4:X e k 

8 i ~ 2 2 ' ( ' 

Depending on the discriminant in (fl9l ), there are three cases to analyze: 

1) If Xn < 4A; /&o> men me discriminant in (fl9l) for each i is non-positive. Recall that 
the stability margin is defined as the absolute value of the real part of the least stable 
eigenvalue, which yields 

5( RPRV ) = |^(s min )| = ^. 

2) If Ai > Akp/bp, then the discriminant in (fl9l) for each i is non-negative, the less stable 
eigenvalue can be written as 

+ _ AA - \/{Xi>b ) 2 - AXeko 2ko 

$t ~ 2 b + y/lft -Ako/Xt 

The least stable eigenvalue is achieved by setting A^ = X N , then have the stability margin 

s^) = msmm)l= f° . 

b o + V b l ~ 4k /X N 

3) Otherwise, if the discriminant in (fT9l ) is negative for small i and positive for large i, then 
the stability margin is given by taking the minimum of the two cases above. This completes 
the proof. ■ 

We are now ready to present the proof of Theorem [Q 

Proof of Theorem [7] We see from Lemma CD that the smallest and largest eigenvalues of matrix 
play important roles in determining the stability margin. To get a lower bound of the 
stability margin, a lower bound for the smallest eigenvalue and an upper bound for the largest 



eigenvalue is needed. Recall that Ai = 2 — 2y/l — e 2 cos#i, Xn = 2 — 2y/l — e 2 cos 9 n, where 

*i E Mwy. 2(ifeiy)> e * E (f^TTf ' ©)• We therefore have X ^ 0, 6 N -+ tt as JV ^ oo, 
and consequently, 

inf Ai = 2-2V1 - e 2 , (20) 

JV 
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sup \ N = 2 + 2Vl-e 2 . (21) 

TV 

To prove the result with RPAV feedback, we consider the following two cases: 

1) Case 1: Ai > bl/4k . According to LemmaQ] the stability margin is given by S'( RPAV ) = b /2. 

2) Case 2: Ai < bl/Ak . From Lemma [U the stability margin is given by 

^(rpav) = ~ V^o - 4fc Ai 



Since Ai > 2 — 2yi — e 2 , we obtain 



^(RPAV) > V_^ _^ i_ (2 2) 

Notice that the above lower bound (|22j) is smaller than &o/2, the value of S^ RPAyS) in case 1. The 
real part sign 3ff(.) in © comes from combining the above two cases. We obtain the first result 
of the theorem. 

To prove the result with RPRV feedback, we consider the following three cases: 

1) Case 1: A^ < 4& /&o- According to Lemma [Q the stability margin is gO^) = 6 Ai/2. 



Moreover, from (l20l) . we have inf/v Ai = 2 — 2y/l — e 2 , therefore the stability margin has 
the lower bound 

2) Case 2: Ai > 4A; /5 2 I . From Lemma [H the stability margin is given by 

Cji(rprv) 



2k Q 



b + ^/b 2 -4k /\ N 



In addition, we have from (f2~Tj) that sup^ A at = 2 + 2vl — e 2 , so the stability margin for 
this case is bounded below as 

g(RPRV) > 

~~ b + ^b 2 -2k /(l + VT^) 
3) Case 3: Otherwise, the stability margin are bounded below by the minimum of the above 
two cases. 

Notice that in the second case, 2k ° > ^. Combining the above three cases, we 

have that 

h 



5 (rpr V ) > min | 6o(1 _ v /r372 ); 

which completes the proof. 
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Fig. 3. Stability margin comparisons between symmetric control and asymmetric control. 



A. Numerical verification for 1 -D vehicular platoon 

In this section, we present numerical verification of the lower bounds of the stability margins 
for both RPAV and RPRV feedbacks with asymmetric control, which are predicted by Theorem [Q 
In addition, the stability margins with symmetric control are also computed to compare with the 
asymmetric case. The stability margins are obtained by numerically evaluating the eigenvalues 
of the state matrix A( RPAV orRPRV ) of © with corresponding controllers. Figure [3] depicts the 
comparisons between the stability margins with symmetric and asymmetric control for the two 
types of feedback: RPAV and RPRV. For both symmetric and asymmetric controls, the nominal 
control gains used are k = 1, b = 0.5, and for asymmetric control, the amount of asymmetry 
is e = 0.1. We can see from Figure [3] that the stability margin of the vehicular platoon with 
asymmetric control is indeed bounded away from uniformly in N, and the predictions Eq. C7]) 
and Eq. © of Theorem Q] are quite accurate. Furthermore, for the same N, the stability margin 
with asymmetric control is much larger than that with symmetric control, especially when N is 
large. 
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IV. Stability margin with D-dimensional lattice information graph 

In this section we analyze a more general scenario than the 1-D platoon of the previous 
sections. We consider a vehicular formation in which the position of each vehicle has dimension 
higher than one, such as a vehicular formation moving in 2-D or 3-D space. We assume the 
dynamics of each of the coordinates of a vehicle's position are decoupled and each coordinate 
can be independently controlled. Under this fully actuated assumption, the closed loop dynamics 
for each coordinate of the position can be independently studied; see ED, EH for examples. The 
information used by a vehicle to compute its control is based on relative measurements with a set 
of neighbors specified in terms of an information graph. The problem formulation is similar to 
the 1-D case in the sense that each vehicle has to maintain constant separation with its neighbors 
in an information graph, except that the information graph now is a D-dimensional lattice. 

Definition 2: An information graph is a graph G = (V, E), where the set of nodes (vehicles) 
V = {1, 2, . . . , N, N + 1, . . . , N + N r } consists of N real vehicles and N r "fictitious" reference 
vehicles. Two nodes i and j are called neighbors if G E, and the set of neighbors of i are 
denoted by jVj. □ 
In this paper we restrict ourselves to D-dimensional lattices as information graphs: 

Definition 3 (D -dimensional lattice): A D-dimensional lattice, specifically a n\ xn 2 x ■ ■ • xno 
lattice, is a graph with nin 2 . . . n D nodes, in which the nodes are placed at the integer coordinate 
points of the D-dimensional Euclidean space and each real vehicle connects to vehicles which 
are exactly one unit away from it. □ 



X-2 



6 5 

) Q- 



^ JL3 

o 0- 



^ 2 I 1 9 

o ch m- 



xi 



Fig. 4. A pictorial representation of a 2-D information graph. The filled node represent the reference vehicles and the solid 
lines represent edges in the information graph. 
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Figure @] depicts an example of 2-D lattice. A D-dimensional lattice is drawn in W D with a 
Cartesian reference frame whose axes are denoted by xi,X2, ■ ■ ■ We also define Nj, (d = 
1, . . . , D) as the number of real vehicles in the Xd direction. Then we have AiA<2 • • ■ Nd = iV 
and riin 2 . . .n D = N + N r . An information graph is said to be square if Ni — N 2 — • ■ ■ — N D . 
Note that the information graph for the vehicular platoon considered in the previous sections is 
a 1-D lattice with N real vehicles (nodes) and N r (= 1) reference vehicle. 

For the ease of exposition, we only consider the case where the reference vehicles are arranged 
on one boundary of the lattice. Without loss of generality, let it be perpendicular to the X\ axis, 
see Figure |4] for an example. This arrangement of reference vehicles simplifies the presentation 
of the results. Arrangements of reference vehicles on other boundaries of the lattice can also be 
considered, which does not significantly change the results; see 11231 . Il24ll . 

Due to its similarity with the 1-D case, we omit the details on desired separations etc, which 
are available in Q. The control laws with RPAV and RPRV feedback, in terms of the errors pi 
are, respectively 

D D 

Ui 

d=l d=l 
D D 

Ui = - k^id +) (pi - pid+) - k^id-)ipi - pid-) 

d=l d=l 
D D 

d=l d=l 

where i d+ (respectively, i d ~) denotes the neighbor of i on the positive (respectively, negative) 
Xd axis. The closed loop dynamics are again represented as x = 

^(RPAVorRPRV)^ where ±Q 

x := \pi,pi, ■ ■ ■ iPn.Pn] ^ K 27V is a vector of the relative positions pi and relative velocities pi. 
The stability margin is defined as before. 

It is shown in that asymmetry in control gains can improve the stability margin with 
RPAV feedback, but the analysis is limited for e — > and the case with RPRV feedback was 
not considered. In this paper, we consider the following homogeneous and asymmetric control 
gains that introduce asymmetry only in the xi axis: 

kufi+) = (1 + e)k , hup-) = (1 - e)k , 
RPAV: (25) 

/c (M d + ) = k , (d > 1), h = b . 



~ k (i,i d +)(Pi - Pi d +) ~ k (hi d ~)(Pi ~ Pi*-) - hp h (23) 
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= (1 + e)ko, = (1 - e)fc , 

RPRV: 6 (M i +) = (1 + e)6o, = (1 - e)6 , (26) 

— = &o, > 1). 

Again, we comment that heterogeneity in control gains has little effect on the scaling trend of 
stability margin, not only for vehicular formation with lattice graphs but also for general graphs 
with bounded degree and weights, please refer to Lemma [2| given in the appendix. 
We first summarize the results in 0, [|23l . 

Proposition 2: Consider a vehicular formation whose information graph is a D-dimensional 
lattice. With the control gains given in (1251) and (1261) respectively. 

1) [Theorem 1 of 0, Theorem 4 of [|23]|] With symmetric control (e = 0), both ,S (RPAv) and 

2) [Theorem 2 of 0] With the control gains given by (1251) . the stability margin with RPAV 



feedback is S (RPAv) = O(^), which hold in the limit e -> and Ni ->• oo. □ 
We next state the main result of this section, which is a corollary of Theorem \T\ It describes 
the stability margin for a vehicular formation with D-dimensional lattice information graph with 
asymmetric control. 

Corollary 1: With the control gains given in (1251) and (|26l) respectively, and < e < 1, the 
stability margin of the vehicular formation with RPAV or RPRV feedback is bounded away from 
0, uniformly in N. Specifically, 

5 (RPAY) >^ V (27) 



S {RPRV) >mmh (l-VT^), 

L 6 J 



(28) 



□ 

Remark 2: From Proposition [2l we see that with the particular arrangement of the reference 
vehicles as mentioned before, the stability margin of the vehicular formation with symmetric 
control only depend on N%, the number of real vehicles along the x\ axis of the information 
graph. For a square information graph, no matter how large its dimension D is, the loss of 
stability margin with increasing number of vehicle N is inevitable, since N\ = N l l D . To make 
the stability margin independent of N with symmetric control, one needs to employ a non- square 
information graph, such that N\ is a constant regardless of the increasing of N . The price one 
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pays is either long range communication and/or increased number of reference vehicles; see 0, 
E3l for more details. In addition, for the RPAV feedback case, with vanishingly small amount 
of asymmetry, the stability margin is improved to 0(1/ Ni), compared to the 0(1/ Nf) trend in 
the symmetric case. 

In contrast, Corollary \T\ shows that with judicious asymmetric control, the stability margin can 
be made independent of the number of vehicles iV in the formation, without using the non-square 
information graph aforementioned. Note that the result we establish in this paper (Corollary [B is 
stronger than that in O, even though the control law is the same. The reason is that the analysis 
in B21 relied on a perturbation technique, which limited its applicability to vanishingly small e. 
In this paper we do not use perturbation techniques, and obtain result for any non-vanishing 
e G (0, 1). In addition, we also consider the RPRV feedback case, while [2J analyzed only RPAV 
feedback. □ 

Proof of Corollary [7] With the control gains specified in (|25T) and (l26l) respectively, it is straight- 
forward - through a bit tedious - to show that the state matrices A^^ and A( RPRV ' can be expressed 
in the following forms, 

= I N ®A 1 + ® A 2 , 
^(rprv) = In ^ A3 + l (d) A ^ (29) 

where Ai, A 2 , A3, A$ are given in (flOl) and has the following form: 

= I Nd ® L^" 1 ) + ® In^.-.n^, 2<d<D, (30) 

where is given in (ITTb and is a matrix of dimension N d x N d , which is given by 

1 -1 



T (d) 



-1 



-1 



-1 



2 -1 
-1 1 



(31) 



The eigenvalues of are given by (see ll2Tll ): 

(id - 1)7T 



2 - 2 cos ■ 



1,2,..., N d . 



(32) 
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For example, for a 2-dimensional information graph shown in Figure HI 



o 


— 1 -I- f 
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-1 + e 


-1 
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-1 - e 


3 + e 





-1 








-1 
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-1 + e 











-1 


-1 -e 


2 + e 



It's straightforward to show that = I 3 <g> + T( 2 ) <g) /2, where T^ 2 ^ is a matrix with 
dimension 3x3. 

From the proof of Lemma [Q we see that the eigenvalues of and ^4( RPRV ) are given by the 

roots of the characteristic equations s 2 + b s + k X^= and s 2 + &oA^s + k X^= respectively, 
where A^is the eigenvalue of L( D \ and £ = (£ 1 , ■ ■ ■ ,£ D ) in which £ d E {1,2, ■ ■ ■ , N d }. We next 
claim that the eigenvalues of D- D > are given by 

D 

\ r =\ h (LU) + J2K(TW). (33) 

d=2 

We prove by induction method. For the case d = 2, = I N , 2 <S)L^ +T^<S)Ini ■ Following (fl6l) 
in the proof of Lemma \T\ the eigenvalues of are given by 

= A £l (L( 1 )) + A, 2 (T( 2 )), 
Now, we assume the general formula for the eigenvalues of L^ D ^ l "> is given by 

D-l 

X iu ... /D _ 1 ^X h (L^) + J2K(T {d) )- (34) 

d=2 

For the case d = D, the matrix has the form given in (|30l) . use (fl6l) again, we have 

^i,...,^ = (J {o-(L (£,_1) + Ai 1J /jv 1 ...jvb_i)} 

= a, 1 .., d _ 1 (l^- 1 )) + a £d (t^)), 

which proves the claim. Now, use (fl4l) and (1321) . the smallest eigenvalue of is equal to Ai, 
the smallest eigenvalue of The result now follows from Lemma [Hand Theorem [T] ■ 
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Fig. 5. Stability margin comparisons between symmetric control and asymmetric control for a vehicular formation with 2-D 
square information graph. 



A. Numerical verification for vehicular formation with D -dimensional information graph 

In this section, we present numerical verification of the theoretical predicted lower bounds 
of stability margin for vehicular formations with D-dimensional lattice information graphs. For 
simplicity, we take 2-D lattices as examples. We assume the information graph is square, i.e. 
Ni = N 2 = y/N. In addition, the stability margins with symmetric control are also computed to 
compare with the asymmetric case. For both symmetric and asymmetric controls, the nominal 
control gains used are k = 1, b = 0.5, and for asymmetric control, the amount of asymmetry 
used is e = 0.1. We observe from Figure |5] that, with asymmetric control, the stability margin of 
the vehicular formation with RPAV or RPRV feedback is indeed uniformly bounded below by the 
prediction Eq. (|27l) and Eq. (1281) respectively. Furthermore the stability margin with asymmetric 
control is much larger than that with symmetric control for the same N. 
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V. Summary 

We studied the stability margin of vehicular formations on lattice graphs with distributed 
control. The control signal at every vehicle depends on the measurements from its neighbors 
in the information graph, which is a D-dimensional lattice. Inspired by the previous works fl}, 
0, we examined the role of asymmetry in the control gains on the closed loop stability margin. 
We showed that with judicious asymmetry in the control gains, the stability margin of the 
vehicular formation can be bounded away from uniformly in N. This eliminates the loss 
of stability margin with increasing iV that is seen with symmetric control. In this paper, the 
analysis of the stability margin avoids the PDE approximation and perturbation method used 
in flTJ , 0. In particular, the latter limited the analyses in those papers to vanishingly small 
amount of asymmetry and resulted a 0(1/ N) scaling trend of stability margin. In addition, the 
control laws examined in JD, required vehicles to have access to the desired velocity of 
the formation. We generalized the results to the case when only relative velocity and relative 
position measurements are available. We showed in this paper that in both cases (i.e., with or 
without absolute velocity feedback), stability margin can be made independent of the size of 
the formation with asymmetric control. The issue of sensitivity to external disturbances with 
asymmetric control is a topic of future research. 
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Appendix 

Lemma 2: Consider a vehicular formation whose information graph is an arbitrary connected 
graph (with 1 grounded node/leader) with bounded degree and weights. With the double integrator 
dynamics £0) and the following heterogeneous and symmetric control gains with RPAV and RPRV 
feedback respectively, 

Ui = - k o w i,j(Pi - Pi)- b oPi, 



(rprv) 



^2 k oWij(pi - pj) - ^2 b oWi,j(Pi - Pj), 



where k ,b are positive constants and w i: j = Wj ti > are the weight on the edge and 
(J, i) respectively. The stability margin decays to zero as N goes to infinity, i.e. 

£ ( RPAV or RPRV) ^ ^^00, 

where N is the number of vehicles in the formation. □ 

Proof of Lemma [2] The state matrices A^ Rmv ^ and A( RPRV ) can be expressed in the following 
forms, 

= I N ® Ar + L® A 2 , 
^(rprv) = In ®a 3 + L®A 4 , 



where Ai,A 2 ,A 3 , A^ are given in (1101) and L is the grounded graph Laplacian of the vehicular 
formation. First, we recall that the Laplacian matrix C of an arbitrary graph G with edge weights 



Wij is defined as 



EfcLi w i,k i = j, (i, k) e E, 



otherwise. 

The grounded graph Laplacian L is obtained by deleting the row and column of C corresponding 
to the leader (recall that we have a leader/reference vehicle). Without loss of generality, let it be 
indexed by 1. The neighbors of the leader is denoted by N\. 

From the proof of Lemma [T] we notice that the stability margin of the formation is determined 
by Ai, the smallest eigenvalue of its grounded graph Laplacian L. 
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We first claim that as long as its grounded graph Laplacian L is symmetric, i.e. Wi j = Wj t i, the 
smallest eigenvalue Ai of L satisfies 

Ai -> as N -> oo. 

This fact can be seen as a generalization of Lemma 3.2 of ll25l . First of all, the grounded graph 
Laplacian L is positive definite 11261 . From the Rayleigh Ritz Theorem E71 . for an arbitrary 
vector x, we have the following inequality 

x T \i(L)x < x T Lx, =>■ Ai < — = — . 

x 1 x 

In particular, we pick the following vector with dimension (N — 1) x 1, x = [1, 1, • • • , 1] T , then 
we have 

x T Lx |7Vi| maxjejvi Wj,i 
N-l ~ N-l 

where \J\fi\ denotes the number of neighbors of the leader. Since the weight Wj t i and are 
bounded, when N goes to infinity, the fact follows. With this, the rest of the proof follows from 
Lemma [Q ■ 
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